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<]} , ABSTRACT. We use free probability techniques for computing spectra and Brown measures of 

f""^ ■ some non hermitian operators in finite von Neumann algebras. Examples include u n + where 

u n and are the generators of Z„ and Z respectively, in the free product Z„ * Z, or elliptic 
elements, of the form S a + iSp where S a and Sp are free semi-circular elements of variance a 
and (3. 
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1. Introduction 



< 

o. 

Recently Haagerup and Larsen QHL99Q have computed the spectrum and the Brown measure 



of i?-diagonal elements in a finite von Neumann algebra, in terms of the distribution of its 
radial part. The purpose of this paper is to apply free probability techniques for computing 
spectra and Brown measures of some non-hermitian, and non-i?-diagonal elements in finite 
von Neumann algebras, which can be written as a free sum of an i?-diagonal element and an 
element with arbitrary ^-distribution. Motivations for this study are twofold, on one hand some 
of these elements appear as transition operators of random walks on groups or semi-groups, see 



. e.g. QdlHRV93aj ], [ |dlHRV93b| ], [ BVZ97Q , here we shall for example treat linear combinations 



of u n and u^, the generators of Z n and Z in Z n * Z and u 2 + v 2 + ?w On the other hand 
random matrix theory has a close connection with free probability (see [ |VDN92P ), but for the 



Q\ ■ moment very little has been done for understanding limit distributions of spectra of non-normal 

matrices in terms of free probability. For example, the empirical distribution on the eigenvalues 
of a random matrix with independent identically distributed complex entries, suitably rescaled, 
converges, with probability one, as its size grows to infinity, to the circular law (the uniform 
distribution on the unit disk), see [ pir84p , QGir97a| ], [ Bai97| ], which is the Brown measure for 



a circular element, in the sense of Voiculescu. It is known that the circular element is the limit 
in * -distribution of the above random matrices, but it is not possible to deduce from this the 
convergence of the empirical distribution on the spectrum (see Lemma 2.1 below). 

Another example that we shall consider in this paper is the free sum of an arbitrary element 
with a circular element. Hopefully, the corresponding Brown measures should represent limit 
of eigenvalue distributions of random matrices of the form A + W where A is a matrix with 
some limit ^-distribution, and W is a matrix with independent entries. In addition to the circular 
element discussed above, this is known to be true for the so-called elliptic element, which can be 
written as S a +iSp and whose Brown measure was first computed in [ Lar99Q by ad-hoc methods. 
It turns out to be treatable by our method as well. The empiciral eigenvalue distribution of its 
matrix model with Gaussian random matrices is computed in [ HP98Q and shown to converge to 
the uniform measure on its spectrum, an ellipse. 

However in this paper we shall stick to the purely free probabilistic aspects of the subject, 
and not touch upon the random matrix problem. We hope to deal with this somewhere else. 
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This paper is organized as follows. In section 2 we recall preliminary facts about Brown 
measures and free probability theory. In section 3 we give a general approach towards the com- 
putation of the Brown measure for the sum of an i?-diagonal element with an arbitrary element. 
We specialize in sections 3 and 4 to the cases where the i?-diagonal element is a Haar uni- 
tary or a circular element, respectively. We close with some final remarks in section 5. The 
pictures of random matrix spectra appearing in various sections of this papers were computed 
with |GNU| |octave| and plotted with |gnuplott the plots of densities of various Brown mea- 



sures, which accompany or replace the rather unwieldy density formulae, were computed by 
|Mathematica| . 
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Institute in Vienna in spring 1999, organized by the Institut fur Funktionalanalysis of the uni- 
versity of Linz and its head J.B. Cooper. The second author was supported by the EU-network 
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2. PRELIMINARIES 

2.1. The Fuglede-Kadison determinant and Brown's spectral measure. Let M be a finite 
von Neumann algebra with faithful tracial state r and denote, for invertible a E Ai, A(a) = 
e *r(iog |a|) j ts F U gi e d e -Kadison determinant (cf. QFK52| ]). Denoting by fx x the spectral measure for 
the self-adjoint element x E M, i.e. the unique probability measure on the real line satisfying 
r(x n ) = J t n dfi(t), we have the following formula for the logarithm of the determinant, which 
serves as a definition of the determinant in the case where a is not invertible 



logA(a) = / logt dfi\ al (t) . 

JR 



The function A (A— a) is a subharmonic function of the complex variable A, and there is a unique 
probability measure ^ on C, with support on the spectrum of a, called the Brown measure of 
a, such that 

log A(A - a) = J log | A - z\ ^ a (dz); 

it is given by 

Pa= ^V 2 logA(A-a) 

where V 2 is the Laplace operator in the complex plane, in the sense of distributions (see 
QBro86| ]). If a is normal, then fx a is just the spectral measure of a. When M. is M n (C), with the 
canonical normalized trace, then fj, a is the empirical distribution on the spectrum of a (counting 
multiplicities). Although the Brown measure of a can be computed from its * -distribution, i.e. 
the collection of all its *-moments r(a £l a £2 • • ■ a 6n ), where a £j is either a or a*, it does not de- 
pend continuously on these ^-moments. Indeed let for example a n be the n x n nilpotent matrix 
with ones on the first upper diagonal and zeros everywhere else, then as n goes to infinity the 
*-moments of a n converge towards those of a Haar unitary (a unitary element u with t{u ti ) = 
for n 7^ 0), whose Brown measure is the Haar measure on the unit circle, whereas the Brown 
measure of a n is 8 for all n. 

Lemma 2.1. Let (a n ; n > 0) be a uniformly bounded sequence whose ^-distributions converge 
towards that of a, and suppose the Brown measure of a n converges weakly towards some mea- 
sure \x, then one has 

(i) J log |A — z\ p(dz) < A(A — a) = J log | A — z\ n a {dz) for all A E C 

(ii) J log |A — z\ fi(dz) = A(A — a) = J log |A — z\ /i a (dz) for all A large enough. 
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Proof. The distribution of |A — a n \ has a support which remains in a fixed compact set, and it 
converges weakly towards that of | A — a\. Part (i) follows from this and the fact that the function 
log is a limit of a decreasing sequence of continuous functions. If A is large enough, then the 
union of the supports of the distributions of the | A — a n \ is away from 0, hence the function log 
is continuous there and (ii) follows from weak convergence. □ 

The outcome of (i) of the Lemma is that the measure /i a is a balayee of measure /i, while we 
get from (ii) the following 

Corollary 2.2. Let U a be the unbounded connected component of the complement of the sup- 
port of n a , then the support of /i is included in C \ U a . 

Proof. The function f log |A — z\ ^ a (dz) is harmonic in U a , while f log |A — z\ /i(dz) is subhar- 
monic there, consequently J log | A — z\ ^ a (dz) — f log |A — z\ n(dz) is a nonnegative superhar- 
monic function on U a . Since this function attains the value by (ii), it is identically by the 
minimum principle, therefore J log |A — z\ fx(dz) is harmonic on U a , and thus the support of fi 
is included in C \ U a . □ 

Conversely, given two measures ji and [i a on C satisfying (i) and (ii), we do not know whether 
there always exists a corresponding sequence (a n ) n >o, fulfilling the hypotheses of Lemma ET 



2.2. R- and ^-transforms. We shall refer to [ |VDN92D , and [ ]Voi98D or PP99H for basic con- 
cepts of free probability theory. Let (A^, r) be as in section 2.1, and let a G M.. The power 
series 



G, 



1 oo 



r(a n ) 



s n=0 S 



can be inverted (for composition of formal power series), in the form 

oo 

K a (z) = - + Y j c n+1 z n = -(1 + R a (z)). 



n=0 



The power series R a is called the R-transform of a (note that this slightly differs from the 
original definition of Voiculescu) and its coefficients are called the free cumulants of a. Let 



iP a (z) = Y^T(a n )z n = l -G a (;) "I 

n=l * ' 



be the generating moment series for a, and assume that the first moment is nonzero, so that 
■0a(O) 7^ 0. Then ^ a has an inverse x a , and the S'-transform of a is defined as 



SJz) 



l + z 



Xa(z) 



Observe that the power series zS a (z) and R a (z) are then inverse of each other (when the mean 
is nonzero). The relevance of these series to free probability is that, if a, b 6 M. are free, then 



Ra+b — R a + Rb 



and 



S a b — S a Si, 



see e.g. QVDN92D . 
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2.3. Calculus of R*-diagonal elements. We use the same notations as in the previous section. 

Definition 2.3. A non-commutative random variable x is called R-diagonal, if x has polar 
decomposition x = uh, where u is a Haar unitary free from the radial part h = \x\. 

Recall that a unitary u E M. is called a Haar unitary if r(u n ) = for all integers n ^ 0. One 
can check that the product of an arbitrary element y with a free Haar unitary is an .R-diagonal 
element. According to QHL99| ], any .R-diagonal element with polar decomposition x = uh has 
the same distribution as a product ah, where h has a symmetric distribution, and its absolute 
value is distributed as h, whereas a is a self-adjoint unitary, free from h, and of zero trace. 
Indeed, one can assume h = a'h, where a' is a symmetry commuting with h and aa' is a Haar 
unitary free from h. Let a, b be two free .R-diagonal elements, then one has equality in *- 
distribution of the pairs (a, b) and (ua, ub) where u is a Haar unitary free with {a, b}, therefore 
a + b has the same * -distribution as u(a + b) which is .R-diagonal, and thus the sum of two 
free .R-diagonal elements is again .R-diagonal. Let f x (z 2 ) = R h (z) be the cumulant series of h, 
which determines the ^-distribution of x, then the power series z(l + z)S x * x (z) and f x (z) are 
inverse of each other. Furthermore if a, b are two free i?-diagonal elements, then one has 

(2-1) fa + b = fa + f b . 

See [|NS97p, [ jNS98| ] and flHL9^ ]. 



2.4. Brown measure of i?-diagonal elements. In QHL99Q the Brown measure of an .R-diagonal 
element is determined as follows. 



Theorem 2.4 ( QHL99| , Thm. 4.4, Prop. 4.6]). Let u, h be *-free random variables in (A4, t), 
with u a Haar unitary and h positive s.t. the distribution \Xh ofh is not a Dirac measure. Then 
the Brown measure fi u h of uh has the following properties. 

{%) fi u h is rotation invariant and its support is the annulus with inner radius Wh^ 1 ]]^ 1 and 
outer radius \\h\\2- 

(ii) The S -transform S^^ of 'h 2 has an analytic continuation to a neighbourhood of "]/x/i({0}) — 1, 
and its derivative S'^ 2 is strictly negative on this interval and its range is «S Mh2 (] A*h({0}) — 
l,0]) = [\\h\\ 2 2 ,\\h-\\ 2 [. 
(Hi) fi uh ({0}) = n h ({0}) and fort e ]p, h ({0}),l] 

Huh \ B I 0, 



(iv) fi u h is the only rotation symmetric probability measure satisfying (Hi), 
(v) Ifh is invertible then a{uh) = supp ji u h, i.e., the annulus discussed above, 
(vi) Ifh is not invertible then a(uh) = B(0, \\hW2). 

The proof involves a formula for the spectral radius of products of free elements. 



Proposition 2.5 ([ |HL99j , Prop. 4.1]). Let a, b be *-free centered elements in M. Then the spec- 



tral radius of ab is 

p(ab) = ||a|| 2 H&H2 

In particular, an .R-diagonal element a = uh can be written as U\U 2 h, with free Haar unitaries 
Mi, u 2 and therefore its spectral radius is p(a) = ||wi||2 1 1 ^2 ^ 1 1 2 = IHh- 



SOME EXAMPLES OF BROWN MEASURES 



5 



3. Adding an ^-diagonal element 

In this section we give a general approach to computing the Brown measure of the sum of a 
random variable with an arbitrary distribution and a free -R-diagonal element. So we let a be an 
arbitrary element, h be self-adjoint and u a Haar unitary, with {a, u, h} forming a free family. 

3.1. The spectrum of a + uh. The spectrum of a + uh is determined as follows. For A G" cr(a), 
A — a — uh is invertible if and only if 1 — uh(X — a)^ 1 is invertible. If h is not invertible, then 
by the result of Haagerup and Larsen on i?-diagonal elements, the latter is the case if and only 
if 

(3.1) \\h(X-a)-%=\\h\\ 2 \\(X-a)- 1 \\ 2 <l; 

if h is invertible, we get the additional possibility that 1 < ||/i _1 || 2 ||A — a||2. In this case we can 
look at (uhy 1 (X -a) - I. 

The case where A G a(a) must be considered individually. Complications arise for such A, 
for which A G a (a), but ||(A — a) _1 || 2 < oo. Otherwise condition ( fOl ) will be satisfied when 
approaching A from the outside of a (a), so that A lies in the closure of the spectrum of a + uh, 
hence in the spectrum. 



3.2. The Brown measure of a + uh. We can assume that u = u\u 2 with u\ and u 2 Haar 
unitaries, where {ui, u 2 , a, h} is a free family, to get 

log A(A — a — uh) = r(log \u*(ui(X — a) — u 2 h)\) 
= r(log |«i(A - a) - u 2 h\) 



J \og\z\ dll ui{ x-a)-u 2 h{.z) 



and this is the Fuglede-Kadison determinant of x\ = ui(X — a) — u 2 h, which is an ^-diagonal 
element whose * -distribution can be computed according to (|2TT|), i.e. f x = f Ul \\- a \ + fu 2 h- 
This in turn will yield the S'-transform of x* x x\, and then by Theorem we can compute 
log A(A — a — uh). 

From the discussion in section ^]3| we have the relation 

(3.2) /ir 1> (o = ^(i+^Q 

In order to be more specific, let us assume that a is self-adjoint, then the computation of the 
distribution of (A — a)* (A — a) is conveniently accomplished by using the Cauchy transform of 
a, namely factoring ( — |A — x\ 2 = (x — x + )(x — xJ) with 

(3.3) x± = ~ ^A + A± ^(X- A) 2 + 4C^ = Re A ± V( ImA ) 2 - C 

and expanding into partial fractions 

1 1 



/° \ ry 2 rf rtr> \ rv* rv* ry 

A Jb\ Ju — \ Ju _j_ Ju Ju — 

we get 

d\i a {x) G a {x + ) - G a (xJ) 



(3.4) C| A _ a | 2 (C) 



( — \X — x\ 2 x + — x. 
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Using the same technique one can compute the 2-norm of the inverse of A — a. Assuming again 
that a is self-adjoint we have that 



||(A-a 



, — 1 II 2 



dfi a (x) 



12 "./ IA-x| 2 

d/J, a (x) 



(3.5) 



(A — x)(X — x) 
1 f ( 1 1 



A — Ay \A — a; A — x 
G a {\) - G a (X) 
A- A 



djJLaix) 



Let us consider the simplest non-trivial random variable, namely a = u 2 = [ ? J ] having 2- 
point spectrum, so that |A — u 2 \ 2 has a Bernoulli distribution. The ^-transform of |A — u 2 \ 2 = 
1 + |A| 2 + (A + X)u 2 is easily computed to be 



R xlxx (z) = (l + \X\ 2 )z+ 1 -[ v /l + l( A -r A) 2 : 2 - I 



and inverting it according to ( |3.2| ) leads to an equation of fourth degree, which apparently is 
unsuitable for further computations. So even this simple case seems to be untractable by this 
method. In fact, so far we have no concrete example where the general method above can be 
carried to the end. We shall develop other methods, in the next two sections, in order to treat 
the cases where the i?-diagonal element is a Haar unitary, and a circular element. 

4. Haar unitary case 
Now a is an element with an arbitrary distribution, free with a Haar unitary u. 

4.1. The spectrum. The spectrum of a + u is determined as follows: one has A £ a {a + u) if 
and only if 1 £ a(u*(X — a)) and since the latter is /^-diagonal, we infer from Theorem |2.4| that 
a necessary and sufficient condition is 



2 



(4.1) \\(X-a)-%i<l<\\(X-a)\\ 
if A £" a (a); otherwise the condition is simply 1 < ||A — a\\ 2 . 

4.2. First approach to the Fuglede-Kadison determinant. We get the following formula for 
the Fuglede-Kadison determinant 

log A(A — a — u) = r(log | A — a — u\) 

= r(log|«*(A-a)-l|) 

log |^ — 1| d/j, u *( X -a)(z)- 



Observe that u*(X — a) is an i?-diagonal element, and we can evaluate the integral as follows. 
The Brown measure of an /^-diagonal element uh is rotationally symmetric with radial distri- 
bution v(dr) and one has 

//"||/l||2 /"27T 
log \z - 1| d/J, u h(z) = / / log \re l9 - 1| d6u(dr) 

JWh- 1 ^ 1 Jo 
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where the inner integral reduces to 



1 f' 2: , 1IM JO r<l, 



— / log|re"-l|d0 
z7T J logr r > 1 



Introduce the radial distribution function 



pr 

F vh {r) = fi uh (B{0,r)) = 2n u{dp) 



which according to Theorem 2.4 is related to the moment generating function ifj h 2 by 



( ^\ ' ' 1 

Kh\r)r 2 



(for || /i" 1 1| 2 1 <r < 1 1 h 1 1 2 ), and by partial integration (note that F( \\ h || 2 ) = 1) 

r\\hh 

r(log \uh — 1|) = / 2-7T log(r) z/(dr) 

JmaxflJI/i- 1 !!" M 



log rF uh (r 



hh _ r m > F uh ( P ) 

iax(l,||h- 1 !|- 1 ) 7 max (l,||/i- 1 ||2 1 ) P 

1 - F uh {p) 



dp 



max(l,||/i- 1 ||2 1 ) P 



dp 



Example 4.1 (2x2 matrix). Let a have the * -distribution of a 2 x 2 matrix, and consider a + u, 
u a Haar unitary. Let p± = p±(X) be the eigenvalues of | A — a\ 2 and let 

...1/1 1 

be its Cauchy transform. Then 

, w = ! G (i)-i = i( _L_ + _L_ 

z \z J 2 \1 — p+z 1 — p_z 

and we get F(r) by solving the equation — t — 1 for t: 

1 1 

2t 



The obvious solution t — 1 is not interesting for us, and dividing it out leads to the other solution 

F(r) 



2(r 2 — p+)(r 2 — 



det A — a 


2 -r 2 r(|A 


— a 


2 ) 


det(r 2 — A — a 


2 ) 
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The logarithm of the Fuglede-Kadison determinant is, for A G a (a + u), 



r(log |A — a — u\ 



-I 



\\X-a\\ 2 i _ 



dr 



+ 



r\\X-a\\ 2 y ( r 

Ji 2 \r 2 - ' r 2 - fi 

](\og\r 2 - +log|r 2 -/i_|) 



dr 



\\X-a\\2 



(log 1 1| A - a\\i - det |A - a\ 2 \ - log |l - 2||A - a\\\ + det |A - a\ 2 \) 



- log 

2 6 



- (log 1 1 -fi+\ + log 1 1 -fi-\) 



It is now convenient to use the representation of the Laplacian in terms of d\ — \ ( d ^ cX — i 
and its adjoint, namely 

V = a(Re A) 2 + <9(ImA) 2 = 4<9 ^ A ' 

Then we have the formulae 



d A ||A-a|| 2 = dAT((A-ar(A-a)) 

= r(A-a*) 
3 A det(A - a) = «9 A ((A - Ax(a))(A - A 2 (a))) 
= 2A - Ai(a) - A 2 (a) 
= 2r(A - a) 



and the density of the Brown measure of a + u is 



(4.3) 

Pa+«(A) = - <9 A <9 A (log 1 1| A - a||* - det | A - a| 2 | - log |l - 2|| A - a|| 2 + det | A - 



7T 



2 f2\\X - af 2 r(A - a*) - 2r(A - a) det (A - a*) 
-Ox 



7i V ||A-a|||-det|A-a| 2 

-2r(A - a*) + 2r(A - a) det (A - a*) 
l-2||A-a|| 2 + det|A-a| 2 

4 ( ||A-a|| 2 - |r(A-a)| 2 | ||A - a||| r(A - a*) - det (A - a*) r(A - a 



7T ^ ||A - a\\j - det |A - o| 2 " (||A - a||* - det |A - a| 2 ) 2 

2|r(A - a)| 2 - 1 |r(A - a*) - r(A - a) det (A - a* 



1 - 2 1| A - a|| 2 + det |A - a| 2 (l _ 2||A - a|| 2 + det |A - a| 
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and in terms of eigenvalues 

Pa+«(A) = -d- x d x ( ~k> 



7T 
7T 



1 



(4.4) 



1 / ~ A*- 



2 

<9a(am 



-(log|l-/z+| + log 1 1 - //-I) 



2 



d\n + + 



7T 



A*+ - 



<9 A (/i+ - A*-) 



- A*- 



+ 



1 / 9x9aA*+ , d\d\[i 



+ 



+ 



<9aA*- 



In particular, if a = [0/3] (Bernoulli distribution) one gets ji± = {|A — a| 2 |, |A — f3\ 2 } and 
consequently the density is 



Pa+«(A) 



|/3 -a 


to 


7T (|A — a 


2 _ 


A-/?| 2 ) 2 



1 

+ 77- 



1 



+ 



1 



A-/?| 2 ) 2 2tt V(l-|A-a| 2 ) 2 (1-|A-/?| 2 ) 



on the spectrum, which is determined by the inequalities 

1 1 



(4.5) 



> 2 



A*+ + A*- > 2 



Specifying further a = 1, /? = — 1, so that a is a symmetry, the spectrum is the region 
bounded by the lemniscate-like curve in the complex plane with the equation 

|A| 2 + 1 = |A 2 -1| 2 

and we get the picture shown in figure [I]. 

This should be compared with the sample fig. |2| of eigenvalues of random 2N x 2N matrices 
of the form X — U2 + Uoo, where U^ is chosen with the Haar measure on U(2N), and U2 = 
VAV*, with V a Haar distributed unitary independent of U^, and A a fixed symmetry of trace 
zero. 

As another example, specify to a = [ [j q ] . As we will see, the spectrum and Brown measure 
are radially symmetric. The eigenvalues of |A — a| 2 are 



(4.6) 



/i± 



t 2 + 2|A| 2 ±tv/t 2 + 4|A| 2 



and hence, substituting this into Q4.5|), we get 



a (a + u) 



X:l4<|A|*< Vl + t + i 



which is a full disk for t > \[2 and an annulus otherwise. For the density we substitute the 
parameters 



I A — a|| 2 



det I A — a\ 



IAI 



into formula ( f43| ) we get the radially symmetric density function 



Pa+uW 



4 



2t 2 



7T V(4|A| 2 + t 2 ) 2 



+ 



(1- 


A 


2\2 




- (1 — 2]A 


2 )t 2 


((1 




A 


2 ) 2 -t 2 ) 2 
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4.3. An alternative expression for the Fuglede-Kadison determinant. In order to treat more 
complicated examples, instead of the integral ( [4.2[ ) it will be more convenient to use a more 
direct formula for the Kadison-Fuglede determinant, which we state as a lemma. 

Lemma 4.2 ([ |HL99| , Proof of Theorem 4.4]). Let uh be an R-diagonal element and define func- 
tions on R + \ {0} by 



f(v) = T((l + vh 2 )- 1 ) 

1 " f(v) 



9\V) 



vf(v) 



Then g(v) is strictly decreasing with g(]0, oo[) =]\\h 1 \\ 2 2 , \\h\\l[ and for every z €]\\h 1 1| 2 2 , \\h\\l[ 
there is a unique v > such that z 2 = g(v). With this v we have 



If 1 

log A(uh-z) = - / log(l + vt) d/i h 2 (t) + - log - 



z 2 



+ vz z 



For our problem of computing log A(A — a — u) = log A(u*(\ — a) — 1) this translates as 
follows. Putting f(v, A) = r((l +v\a — A| 2 ) -1 ) and denoting v(X) the unique positive solution 
of the equation (1 + v)f(v, A) = 1, then 



logA(A -a-u) = log A(>* (A -a) - I) 



- r(log(l + v\a - A| 2 )) - - log(l + v). 



Note that this approach cannot be used in the general setting of section as it does not tell 
how to evaluate the Kadison-Fuglede determinant at z = 0. 

For the rest of this section we shall assume that a is normal with spectral measure fj, a , so that 
we can write 



(4.7) f(v,X) 



and again with (1 + v)f(v, X) = 1, 



dfi a {t) 



1 +v A- t\ 



log A(a + u — X) — - J log(l + v\X - t\ 2 ) dfjt a (t) - - log(l + 
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For the density of the Brown measure we obtain 

2 

p(X) = - dxd\\ogA(a + u - A) 

TV A \J l+ V \X-t\ 2 PW 1+V 

( 2 ^ 

dxv (/ i + v~ix - t\ 2 Mt) - TT^) + "/ i + !|A-tp ^ (t) 



-0* 

7T 



V 3 / 



1 n f 1 f |A-t| 2 i / \ 

a / t — 77- — nr — Tja w) 

7T J X — t 1 + V \X — t\ Z 



1 f 1 



7T / X-t 



\x-t 


2 d~ x 


v + v(X- t) 


(l + v\X-t 


2y 



dfi(t) 



I 0*7 (l + v\\-t\*)* dti{t)+ J (l + v \l-t\>)> Mt) 



now by implicit differentiation 



d x v 



1 = (1 + *)/(«, A) 

= d- x v f(v, A) + (1 + A) d~ x v + d-J(v, A)) 



f+(l + v)d v f 



<KI(r.\) = - i * V( \ %^ d»{t) 



{1 + v 


X- 
-t 


-t\ 2 ) 2 
2 


(1 + v 


X - 


- 1 


2\2 



d v f(v,X) = - / — — - — -r—dfl(t) 



and thus 



(4 - 8) pW = = ( ,if ( ,,» + ('iVw(M)i |Sa/( "' a)|2 + v,(v ' A) + 1,2 s " /( "' A) 



TV 



We will apply this in three situations here. First consider a finite dimensional normal operator 
a, like e.g. a = u n , the generator of the von Neumann algebra of Z n , then the integrals become 
finite sums and can be evaluated numerically. As an example see fig. |3], which should again be 
compared to the corresponding samples of spectra of random matrices in fig. f|. There U3 is a 
fixed 150 x 150 permutation matrix with the same spectral distribution as -u 3 and Uoo is again a 
150 x 150 standard unitary random matrix. 

Secondly, assume that a is self-adjoint. Then we can factorize the denominator in the integral 
QE7t ) as 1 + v \X - t\ 2 = vt 2 - v(X + A) t + 1 + v\X\ 2 = v(t- z )(t - z ) where 



z = Re A + -^/v 2 (ImX) 2 + v. 
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From this we can express f(v, A) and therefore p(X) in terms of the Cauchy transform G(C) of 
a as follows. 

dfi(t) 



v(t- z )(t - Zo) 

1 f 1 ( 1 



V J Zq- Zq \t- Zq t- Z 

Im G(zq) 



y/v 2 (Im\) 2 + v 

As an example consider a = w 2 + t> 2 , where m 2 and w 2 are the generators of two free copies 
of Z 2 . Then a is self-adjoint and distributed according to the arcsine law (or Kesten measure) 
and has Cauchy transform G(Q = — , 1 1 . A picture of the density of the Brown measure of 

u 2 + v 2 + u is presented in fig. q[ 



0.5 



-0.5 




Figure 5. Density of n U2 



Finally, let us consider the free sum of an arbitrary unitary v and a Haar unitary u. Let 
dfi(9) be the spectral measure of v on the unit circle. For the evaluation of the integral Q4.7| ) we 
factorize the denominator again, this time writing 



f(v,X) 



l + v\X-e ie \ 2 

dfi(9) 



l + v{\\\ 2 + 1) - v {Xe- i0 + Xe ie ) 

„i6 



vXe™ - (1 + v{\X\ 2 + l))e ie + vX 

„i8 



dfi(9) 



1 



(e ie - z + )(e ie - z-) 



where 

z± = i (l + -^(1 A| 2 + 1) ± v/(l + ^(|A| 2 + l) 2 ) (1 + ^(|A| 2 - l) 2 )) . 
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Note that = \j\ and \z + \ > \z-\, and thus \z + \ > 1 > 



1 [ 1 f z+ z_ 



vX J z + z^ \{z+ — e %e ) (z- — e te ) 
z + G(z + ) — z_G(z_) 
vX(z + — Z-) 

z + G(z + ) — z-G(z-) 

" TP + ^(i A i 2 + W) (i + v{\x\^rf 

For the determination of the spectrum ( ff~T| ) we need 



dn(6) 



ll(A-«) 



-1||2 



2 / iA-e ie l 2 



|A| 2 - 1 7 V^-e^ A-e-* 9 
AC7(A) + AC7(A) - 1 
~ |Ap^I 

As an example let us consider for q E [—1,1] the unitary u q with Poisson distribution, i.e. whose 
moments are r(w") = gl n L For q = this is the Haar distribution, while for q — 1 it is the Dirac 
measure at 1. By Fourier transform, the density of the spectral measure is 

d N {9) 1 1 _ q2 

The Cauchy transform is 

G q (0 = 

'.C-9 

and from this we get the other relevant functions 

f \M 2 -q 2 




(l-|Ap)|A-<rT 

qz + -q~ 1 Z- 

f{v,X) 



X < 1 



(z + - q)(z- - q- 1 ) v/(l + V (|A| 2 + l) 2 ) (1 + v(|A| 2 - l) 2 ) 
Substituting this into (|4.8|), we get pictures like fig. EL where q = 0.7. 



5. Adding a circular element 

A standard circular element has the * -distribution of C = S\ + 1S2 where Si, S% are free 
standard semi-circular elements, i.e., self-adjoints whose distribution is the semi-circle law 
7^:\/4 — x 2 dx on [—2, +2]. Its polar decomposition is C = uh with u a Haar unitary free 
with h (hence C is i?-diagonal), and h has the quarter circular distribution — x 2 dx on 

[0, Vo]. The symmetrized h in Haagerup-Larsen's decomposition C = ah has a semi-circular 
distribution of variance 2. In this section we consider the Brown measure of X t = X + C t , 
where X has arbitrary * -distribution, it is free with C t and C t is a circular element of variance 

t, i.e. C t = \j\C where C is a standard circular element. It will be convenient to assume that 
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2 







-1 



-2 




-2-1 12 



Figure 6. Density of n Uq+Uoo at q = 0.7 



the C t form a circular process, i.e., for each s < t, C t — C s is *-free with C s . We shall use a 
heat equation like approach, by differentiating in t. One has 

log A(A -X t ) = \ log A(|A - X t \ 2 ) = \ limlog A(|A - X t \ 2 + e 2 ). 

Let us denote H tj£ = \ A — X t \ 2 + e 2 and compute the derivative log A(H t>£ ). To this end let 
dt be small, dC t = C t+dt - C t (so that T(dC;dC t ) = dt). Then 

Ht+dt,s — | A — X t+dt \ 2 + e 2 

= \\-X t -dC t \ 2 + e 2 

= \\- X t \ 2 - (A - X t )*dC7 t - dc;(\ - X t ) + \dC t \ 2 + e 2 

= fr t)C - (A - x t ydc t - dc;(\ - x t ) + dc;dc t 

= H t , e [I - h- £ \(\ - X t )*dC t + dc;(\ - X t ) - dCfdCt)] 

and hence 

log A(H t+dti£ ) = log A(H t , e ) + log A(l - ^((A - X t )*dC t + dC;(X - X t ) - dC* t dC t )) 
= log A(H tte ) + r(log |1 - ^((A - XtfdCt + dC* t {X - X t ) - dC*dC t )\). 
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Now observe that 

r(log |1 + a^Jdi + bdt\) = Jr(log |1 + aVdi + bdt\ 2 ) 



— t (log(l + (a + a*)Vdt + (b + + a* a dt + £>((dt) 3/2 )) 

— t I (a + a )v / rft + (6 + b* + a*a)dt - ^(a + a*) 2 d^) + £>((dt) 3/2 ) 

~t ((a + a*)^ + (b + 6* - Q2 ^ a * 2 )d^ + 0{{dtf 2 ) 



In our situation we have a = -H^ ((A - X t )*^ t + ^i*(A - X t f) and b = H^ 1 ^^, so 
that t(o) = and r(b) = r(if t ~ 1 ) by freeness of dC t and {// t>£ , A — X t }. Further we have 

and using the formula r(a 1 fe 1 a 2 6 2 ) = T(ai)r(a 2 )r(&i& 2 ) if {ai,a 2 } is free from {&i,6 2 } and 
r(bi) = r(6 2 ) = 0, we see that only the last term is nonzero and equal to 

= 2r((A - X t ) H-l (A - X t f) T (H£) 
= 2r(H- e 1 (H t>£ -e 2 )r(H- £ 1 ) 
= 2r(H^) + 2e 2 r(H^) 2 

so that 

l0gA( ^ + "' £ ^ l0gA( ^ £) = \ (2r(H-J) - 2t(H£) + 2e 2 r(H^ + 0{{dtf' 2 ) 

= e 2 r(H- £ r + 0{{dt) 112 ) 

hence 

|logA(^, £ ) = £ 2 r(^ 1 ) 2 



and 



log A^" 1 ) = log m ,e) + f t{H~1? ds. 

Jo 



Let a\ yS be a self-adjoint element with symmetric distribution, whose absolute value is dis- 
tributed as | A — X s \. Now note that by the Stieltjes inversion formula 

^(//- £ 1 ) = r( £ (|A-X s | 2 + £ 2 )- 1 ) 

= -r (lm[(ie - a AiS ) -1 ]) 

dn axs (x) 



,0 7T 



dx 



x=0 

i.e., the density at of the distribution of a\ iS . Now we need the following 
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Lemma 5.1. Let a be a self adjoint symmetrically distributed element, free with S and C, where 
S and C are a semi-circular and a circular element of same variance respectively, then \a + S\ 
and | \a\ + C\ have the same distribution. 

Proof. Let b be a symmetry free with {a, S,C}, then by QHL99| , Prop. 4.2] ba and bS are *-free, 
thus | a + S\ = | ba + bS | is distributed as \ba + C\. Now using the fact that multiplying with a 
free Haar unitary u does not change the * -distribution of C, we can replace the latter according 
to C = u*C, and get the following equalities of ^-distributions 

\ba + C\ = \ba + u*C\ = \uba + C\ = \u\a\ + C\ = \ \a\ + C\ 

□ 



Using the lemma we get 

| A — X 8 \ = | A — X — Cg\ = \a\ + S s \ 
where a x is the symmetrization of | A — X Q \ , free with the semicircular S s and therefore 



\a\ + S s 



It follows from Corollary 3 of QBia97| , p. 711] that the distribution of a\, s has a density at 
which is p s (0) = — , with 



(5.1) 



v(s) = inf < v > : 



dfi\x-x \(x) I 
X 2 + V 2 ~ s 



If A ^ a(X ), then by e.g. [ Bia97 ] 



(JJ-X\^ l^ a \,s ( ) 

er{H ) < sup — -j 

X6R dx 
1 

< 



7T\/ S 



furthermore for s small enough, A ^ cr(X s ) and r(|A — X S A| 2 ) is bounded above, hence 
£t{H~1) also, therefore we can apply the dominated convergence theorem and we get 

log A(A -X t ) = \ limlogA(tf 0)£ ) + \ f e 2 r{H^fds 

2 e^O 2 J 



(5.2) 



lo g A(\-X ) + ±£ 



v(sf 



logA(A-X ) + 



ds 



ds 



where t x = mi{t : v(t) > 0} = (/ 
Brown measure is 



^fJA-Xo|(f)\-l 



. So whenever A ^ cr(X ), the density of the 



1 P t v(s) 2 



7T 



TV 



d\v{s) 



ds 



v(tx) 2 



t 2 

A 



d\t\ 



and the second summand will be zero if is continuous at t\. 
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Example 5.2 (2x2 matrix). Let X = a be as in example 4.1, and consider X t = a + C t . Let 
again fi± be the eigenvalues of (A — a)* (A — a), then the relevant parameters are 



|| A — a 



,2 



2 



IKA-a)- 1 !!!- 1 ^ 1 ■ ^- l|A - a|1 ^ 



2 \fJ,+ fJ,-J det | A — a| 2 



djx\x-a\{x) \ 1 _ ._ lM _ 2 det|A-a| 



2 



The function u (s) 2 is the solution of the quadratic equation 



11/1 1 

+ 



s 2 \fi + + v 2 fi- + v 2 
fJ,+ + /i- + 2v 2 



2 )U, + jU_ + (fM + + yU_)v 2 + V A 

||A - a||| + v 2 
det|A-a| 2 + 2||A-a||^ 2 + i; 4 



which is explicitly 

v(s) 2 = ±(s - 2||A - al^ ± - 2 || A - a|| 2 ) 2 - 4(det |A - a| 2 - s || A - a|| 2 )) 

= l( s ~ 2 II A - a \\l ± ^/s 2 + 4(||A-a||f-det|A-a| 2 )^ 

Now we have to choose the right branch of the square root. To this end, let us compute the 
spectrum of X t : Assume A £ cr(a), then A G a(a + C t ) if and only if 1 — C t (A — a) -1 is not 
invertible. Now C t (X — a)^ 1 is i?-diagonal and not invertible, so by Theorem ^4|(0), 1 is in 
its spectrum if and only if its spectral radius is at least 1 and using Proposition [23] we get the 
inequality 

l<p{C t (X-a)- 1 ) = \\C t \\ 2 \\{X-a)-% 

in other words, 

det|A-a| 2 < t\\\-a\\ 2 2 

and hence for s < t, det |A — a,| 2 — s ||A — «||1 < 0, only the "+" branch gives a nonnegative 
solution. Consequently 

logA(A-A t )-logA(A-A ) 

1 f* 1 ||A-a|| 2 t v /s 2 + 4(||A-a|| 2 1 -det|A-a| 2 ) 



2 J tx 2s s 2 2s 2 



1 -\ogs+ l|A 2g a|1 ^ + \ log ( 8 + \/s 2 + 4(11 A - a\\i - det | A - a\ 2 ' 



{ u / s 2 + 4(||A-a||4-det|A-a| 2 ' 



s=t\ 



Now observe that 

n\\ \ 11/1 1 i I \ 
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and hence, denoting 

(5.3) R{\) = 4(|| A - a\\i - det |A - a\ 2 ) = - fi_) 2 

we get 

logA(A-X t )-logA(A-X ) 

= ilogt + + ^log (t+y/W+m) - ^V* 2 + ^(A) 



1 det |A - a| 2 ||A-o||^ 1 det |A — a| 2 + 2||A — — det |A — a| 2 
4 ° g ||A-a|| 2 ~ 2det|A-a| 2 ~ 4 ° g ||A-a 112 



2 



||A-a||| 2||A-a||t-det|A-a| 



4det|A-a| 2 ||A-o| 



1, , HA-a 
-logt + 



4 ° 2t 
-ilogdet|A-a| 2 -^log2-i 



+ i log (t + v^ 2 + i2(A)) - i Vt 2 + i?(A) 



and finally the density is (note that d\d\ log det |A — a\ 2 = and <9 A <9 A || A — o||| = 1) 
Po+ c t (A) = ^d A logA(A-X t ) 

4^^( log ( t+ v™)-^±M) 

11./' 1 l\ d\R(\) 

— + — C\ 



2tt U+ v /t 2 + J R(A) t 2^/WTRjX) 



|<9 A i?(A)| 2 / , a fl(A)^ A fi(A)-|a A fl(A)| s 

xv^T^ + (* - ^ + ^ (A) J ffi^ 



vrt 4vrf y 2R(X)y/t 2 + R(X) V ) R(Xf 

Again we can specify to a = [§ q], and we get the spectrum 

a(a + C7 t ) = {A: |A 2 -l| 2 <t(|A| 2 + l)} 

Note that for t — 1 this is the same as cr(-u 2 + u) from example 4.1. However this time the 
density is a function of the real part alone, namely substituting y,± = | A ± 1| 2 into (|5.3D , we get 
R(X) = 4(A + A) 2 and consequently the density depends only on the real part 

Pa+cAx + iy) = — + 



7rt 8irx 2 VVt 2 + 16x 2 

The situation for the nilpotent 2x2 matrix a = [q J] is as follows. We have computed the 
eigenvalues of |A — a\ 2 in (|4.6|), and thus 

a(a + C7 t ) = {A:2|A| 4 <t(l + 2|A| 2 )} 



which is the disk with radius y y ^ + ~ + |. This is the same as cr(a + \^tu), but with the 
possible hole removed. Furthermore we get R(X) = — /U-) 2 = 1 + 4|A| 2 and the density 
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function is again rotationally symmetric: 

1 (1+4|A|V* 2 + 1 + 4|AP (1 + 4|A| 2 ) 2 



7Tt 



Example 5.3 (Elliptic law). An interesting example is given by the so-called elliptic random 
variable S a + iSp, where S a and Sp are free semicircular variables of variances a and (3. Note 
that for a = (3 this is a circular variable C 2a - The Brown measure has been computed by 
Haagerup (unpublished) by another method. The name elliptic stems from the shape of its 
spectrum, which is an ellipse. This can be seen as follows. Assuming that a > (3 let 7 = a — f3, 
then for A £ cr(S 7 ) = [—2-fy, 2^7] we have A e cr(S' 7 + C 2 p) if and only if 1 - C 2/3 (A - 5' 7 )~ 1 
is not invertible. From Theorem ^4] we infer that the spectrum of C 2 p(X — 5' 7 ) _1 is the disk 
centered at zero with radius ||C 2( g(A — S , 7 ) _1 || 2 , so that we get 

v(S y + C 2 p) = {A : 1 < 2(3 ||(A - SV)" 1 !! 2 } 



We use formula ( |3~3[ ) for the Cauchy transform Gs 7 (C) = ~ ^ 7 ^ t0 g et 



and hence the spectrum is 

y/A 2 - 4 7 - y/A 2 - 4 7 7+^ = « 

A — A ~ (3 (3 

Now consider the Zhukowski transformation / : £ 1— > | + 7^, which maps the circles {^- : < 
9 < 2ty} to the ellipses { + t) cos + z ( j — t) sin 6 : < < 27r} and hence the open disk 
{£ : |£| < -±=} bijectively onto C \ [—2^/7, 2^/7]. Note that the excluded interval is exactly 

the spectrum of S T So assume that A = /(£) with |£| < is not in the spectrum of S 7 , then 
observe that 

A 2 -47 = i + 2 7 + 7 2 e 2 -47= Q - 7 £ 
and hence A G cr(S' 7 + C 2/ g) if and only if 

a |-|-T^ + 7e_ l+ 7 |ep 



/3 - |-|+7e-7e i-7iei 2 ' 

This inequality reduces to 

a + p 

thus 

a(S 7 + 5^) \ [-2^7, 2^ = : — < |e| < 

and taking the closure of this set we obtain o r (5' 7 + C 2 p) as the interior of the ellipse 

(5.5) ( 2a cosfl+ 213 -.ismO :0 < 6 <2n 
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Now let us turn to the Brown measure. As already noted, the method from section will not 
work on a = S 7 . Indeed the /^-transform of | A — 5* 7 | 2 can be computed from the inverse of 



r , , 1 I, a/4-47 - V^- 4 7 



27 



where x± are as in (|3.3|). Let A = £+ir}, then we can rewrite x± = £± \J (, ~ V 2 an d abbreviating 



V — vC — r ? 2 ' solve the equation G*|a— (C) — z f° r 2/» which gives 

e 2 i 



(1 - 2 7 ^) 2 + z(l - jz) 



It follows that .fir(z) = y 2 + 7] 2 and 



fi M7 p(.)=^(z)-l 



72; £ 2 2 



+ rfz] 



' (1 - 2-fz) 2 

for real A this has been used in [ |HKNY99| ] to characterize the semicircular distributions. In 



order to get the determining series f u \\s \ according to Q3.2D one has to solve a fourth order 
equation, which is not suitable for further computations. So we have to use formula (^2|), for 
which we need v(s) from ( p.l| ) first. We have done most of the work already, since J - ' 

-C|A-5 7 |2(-v 2 ),thus 



\\-x\ 2 +v 2 



v(sf 



-K\ 



|A-5 7 P 



2„2 



and 



(A + A) 5 



's + 27)2 s + 7 
1 (A - A) 2 



+ 77' 



4(s + 27) 2 s + 7 



4s 2 



and the density becomes, with 



P5 Q +i5 /3 (A) = / — = — ds 



7T 



'A 



1 Z' 2 / 9 

-9s 



7T 



2(A + A) 2(A-A) 



4(s + 2 7 )< 



4s 2 



ds 



J_ a _ / A + A _ A- A 
2tt a Vs + 27 s 
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(5.6) 



1/1 I 

+ 



4tt \a /V 2tt A V^a + 27 



1 / A + A _ A^A 



Now t A = || (A — S^) -1 H2 2 has been computed above in ( |5.4|) , and denoting a; = ^A 2 — 47, it 
is 

-1 

I ui — UJ 
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and we claim now that the second summand in ( |5.6| ) is zero. For this note that u 2 — uj 2 = A 2 — A 2 
and hence 

1 n fA + A A-AA 1 . -. / A - A A /x -^foo-u 1 

' 9x((A + A)-(w + £j)-(<d-£j) + (A-A)) 



47T7 

1 <%(2A-2u;) 



47T7 



0; 



Thus we get that the density is constant + jA on the interior of the ellipse Q5.5D 

The elliptic law appears in the random matrix literature in QGir97b| ] . 



6. Other examples 

There are some other examples that can be done by ad-hoc methods. 

Example 6.1. Consider two freely independent symmetries u 2 and v 2 of trace zero, for example 
the generators of the left regular representation of Z 2 * Z 2 . Here we compute the Brown measure 
of T = au 2 + (3v 2 . To get its spectrum, look at its square 

(cra 2 + (3v 2 ) 2 = a 2 + f3 2 + a(3(u 2 v 2 + v 2 u 2 ) 

Since u 2 v 2 = (v 2 u 2 )* is a Haar unitary, we see that T 2 is a normal element with spectrum 
a{T 2 ) = a 2 + 1 + a(3[-2, 2]. Since T and — T have the same distribution, it follows that 

a{au 2 + (3v 2 ) = a 2 + [3 2 + apt : t G [-2, 2] 

The Brown measure can be deduced by the same symmetry considerations, but for the sake 
of simplicity let us consider the special case a = 1, (3 = i only. Here the spectrum is the 
union of the complex intervals [—1 — i, 1 + i] and [—1 + i, 1 — i). The Brown measure of 
(u 2 + iv 2 ) 2 = i(u 2 v 2 + v 2 u 2 ) is the arcsine law (we are taking the real part of a Haar unitary) 

, . . dt 
dv(t) = -j=== 

7TV4 - t 2 

on the imaginary axis. By symmetry considerations we must have the same measure on each of 
the four "legs" of the spectrum, call it n , which must satisfy 

1 r 2 dt 
f(t 2 )dMt) = ~ / f{t) f -— 2 
£ J 7TV4 - t z 

^ u 

f{u 2 ) — -^=^= du 

lo 7rv4 — u A 

and it follows that the density of the Brown measure is 

da ( — ■=- 1] = du (\t\) = — = dt 

Example 6.2. Other examples that are perhaps attackable arise from the following matrix mod- 
els. Consider U 2 + A, where U 2 E U (2N) is a unitary matrix s.t. U 2 = U 2 and tr U 2 = 0, while 
A is an arbitrary 2N x 2A^ matrix. The spectrum of U 2 + A can be bounded as follows. As- 
sume x is a unit eigenvector of U 2 + A with eigenvalue A, then it can be decomposed along the 
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spectral projections of U%: x = x + + x_ so that U 2 x = x + — X-. By assumption we also have 

(U2 + A)(x + + X-) = X(x + + X-), and thus 

x + = ~(1 + A — A)x X- = -(1 — A — 

now by orthogonality (x + , xJ) = we get 

= ((1 + A — (1 - A + A)x) 
= (l + \)(l-X)\\x\\ 2 + (l + X)(x,Ax)-(l-X)(Ax,x) - \\Ax\\ 2 
= (1 + A - A - |A| 2 )||x|| 2 + (A + l)(Ax,x) + (A - l)(Ar, x) - \\Axf. 

Separating real and imaginary part results in two equations 

1- |A| 2 - \\Ax\\ 2 + X(Ax,x) + X(Ax,x) = 
A — A + (Ax, x) — (Ax, x) = 

Let us no consider two specific cases. 

A is unitary: In this case \\Ax\\ = 1 and p = (Ax, x) satisfies the following equations 

-|A| 2 + Ap + Ap = 
A — A = p — p 

or in other words 

|A-p| 2 = |/f 
Im A = Im p 

thus A — p is real and we have 

X-p = ±\p\ 

i.e., 

X e {p±\p\ : p = (Ax, x) G co a (A)} 
A = iB is purely imaginary: Here we assume A + A* = and the equations are 

1- |A| 2 - \\Bx\\ 2 + i(X- X)(Bx,x) =0 
A — A = 2i (Bx, x) 

Hence 

Im A = (Bx, x) 
(Re A) 2 = 1 - ||£x|| 2 + (Bx,x) 2 

If one puts A = UU3U*, where U3 is an 6iV x 6iV model of the generator of Z 3 , and U is a 
random unitary 6A^ x 6iV matrix, then possible eigenvalues are enclosed by the region shown in 
figure [7[ And indeed, samples of small numeric random unitary matrices U2 + UU3U* have an 
eigenvalue density as shown in figure ^J, while in bigger dimensions the eigenvalues concentrate, 
cf. figure |9|. We were able to compute the spectrum of the free sum u 2 + u 3 recently and will 
investigate this topic further in future work. 
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Figure 7. Possible spectra of random U 2 + U 3 
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Figure 8. 5000 samples of eigenvalues of 6 x 6 random matrices U 2 + U 3 
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Figure 9. 200 samples of eigenvalues of 150 x 150 random matrices U 2 + U 3 
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